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ABSTRACT 

We present spatially and spectrally resolved Bry emission around the planet-hosting, 
transitional Herbig Ae/Be star HD 100546. Aiming to gain insight into the physical 
origin of the line in possible relation to accretion processes, we carried out Bry spectro- 
interferometry using AMBER/VLTI from three different baselines achieving spatial 
and spectral resolutions of 2 - 4 mas and 12000. The Bry visibility is larger than 
that of the continuum for all baselines. Differential phases reveal a shift between the 
photocentre of the Bry line -displaced 0.6 mas (0.06 au at 100 pc) NE from the 
star- and that of the AT-band continuum emission -displaced ^ 0.3 mas NE from the 
star. The photocentres of the redshifted and blueshifted components of the Bry line 
are located NW and SE from the photocentre of the peak line emission, respectively. 
Moreover, the photocentre of the fastest velocity bins within the spectral line tends 
to be closer to that of the peak emission than the photocentre of the slowest velocity 
bins. Our results are consistent with a Bry emitting region inside the dust inner rim 
(< 0.25 au) and extending very close to the central star, with a Keplerian, disc-like 
structure rotating counter-clockwise, and most probably flared 25°). Even though 
the main contribution to the Bry line does not come from gas magnetically channelled 
on to the star, accretion on to HD 100546 could be magnetospheric, implying a mass 
accretion rate of a few 10“^ Mq yr“^. This value indicates that the observed gas has 
to be replenished on time-scales of a few months to years, perhaps by planet-induced 
flows from the outer to the inner disc as has been reported for similar systems. 

Key words: stars: individual: HD 100546 — Protoplanetary discs - Line: formation 
- Accretion, accretion discs - Techniques: interferometric 


1 INTRODUCTION 


The Bry emission line at 2.166167 /rm is commonly 
used as a quantitative accretion indicator in pre-main se- 
queii c e (PMS) stars llM uzCTolle et al.l Il998l: I Calvetetal] 


2004; iDonehew fc BrittairJ ^lll : Mendigutfa et al. 20 111 . 
20130 . The importance of Bry as an accretion tracer re¬ 
lies on its ubiquity in young stellar objects, and on the 


* Based on observations made with ESO Telescopes at the La 
Silla Paranal Observatory under programme ID 093.C-0339 
f E-mail: I .MendigutiaOleeds. ac .uk 


fact that the near-IR is much less affected by extinc¬ 
tion than direct accretion-related UV signatures. Indeed, 
line modelling is able to reproduce observed Bry profiles 
mainly in terms of accr e tion/ w inds (iMuzerolle et al. l 200lJ; 
iKurosawa fc Romanoval l2013l : iTambovtseva et al.l 1201^ . 
Due to their higher brightness and sizes when compared 
with classical T Tauris (CTTs), the inner discs of nearby 
Herbig Ae/Be (HAeBe) stars can be spatially resolved 
through interferometry, providing direct observational con¬ 
straints on the nature of the Bry emission. Near-IR spectro- 
interferometry indicates that for most HAeBes the Bry 
emitting region is more compact than the continuum emis- 
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sion arising from the dusty inner disc. For some of these 
cases, the line has been d i rectly associated w i th acc retion in¬ 
flows (jKraus et ahl l2008l : lEisner et al.l l201fll . l2014h , but not 
in others where the line emitting region is more extended 
jMalbet^^^ 2007; iTatulli et al.ll200^ : iKraus et al.ll200^ : 
IWeigelt et al.l 2011 ). HAeBes have an additional nterest 


given that they represent a class of objects where the accre¬ 
tion paradigm could change from magnetically channelled 
such as in CTTs (magnetospheric accretion; MA), to a di¬ 
rect di sc-to star accretion through a boundary layer (BL ; 
see e.g. Blondel fc Tiin A Diidlioo^ : iMendieutfa et al]l201ll : 


ICaulev fc Johns- Krullll2014ll. 

HD 100546 is a nearby (d ~ 100 pc; Ivan den Ancker et al.l 
1 19981 : Ivan LeeuwenI l2007l l transitional HAeBe (B9) star 
with a complex circumstellar environment. It shows an in¬ 
ner dust disc from ~ 0.2 to 4 au, a gap from ~ 4 to 
13 au, and an outer disc extending several hun dreds of 
au (see e.g. iBenistv et J]|2010l: iTatulIi et al.ll20lll. TATll 
hereafter). Molecular gas traced by C O, OH and CH~*~ 
accumulates mainly i n the outer disc (|yan_deLj?las_eLal 


200£; ITM et ^ 1201 ll : iLiskowskv et al.l l2012l : iBrittain et al 


20130 . but there is evidence of atomic, low density gas 


traced by the [O I] 6300 line com ing from the same spa¬ 
tial s cale as the inner dust disc (lAcke fc van den Anckeil 
l2006h . The presence of variable emission lines in the opti- 
cal/UV has been interpret ed as the signature of MA/winds 
operating in HD 100546 llVieira et al'l Il 999l: iDeleuil et al] 


|2004 [Ouimaraes et al.ll2o'o^ IPogodin et al. f l2012ll . The fact 

that there is compelling evidence that HD 100546 hosts 
at least two sub-stellar/planetary objects makes this sys¬ 
tem t he subject of intense stu d y (see e.g. BouwmaJi_eLalJ 
200^;_ Anke_&vanden_^ncke^ 20061: Quanze£_al|_ 201^ 


Mulders et al.ll2013l : IBrittain et al.ll201E , 20141 : IWalsh et akl 

20141 ). In particular , a protoplanet is l ocated in the outer 


UI.J.CI/X U iO ± W V-/Ct U V-l. i.AJ. UAAV./ 

disc, at ~ 50 au iouanz et al. I I 2 OI 5 I . I 2 OI 3 I : ICurrie et al.l 
|2014^. and another candidate in the disc gap at ~ 10 au 
I Brittain et al. I [ 2 ^. 

Aiming to gain insight into the properties of atomic gas in 
the innermost regions of the disc, and their possible rela¬ 
tions with accretion processes in HAeBe stars, we present 
the first high spatial resolution spectro-interferometric data 
of HD 100546, centred on the Bry line. Along this pa¬ 
per, the most recent ste llar parameters deriv ed from a X- 
Shooter/VLT spectrum (iFairlamb et al.ll2015l l will be con¬ 
sidered: T* = 9750 ± 250 K, M. = 1.9 ± 0.4 Mq, R. = 1.4 
± 0.5 R(T). These are consistent with previous determina - 
tions (Ivan den Ancker et al.lll998l : iGuimaraes et al.ll2006ll . 
Section [2] describes the observational set-up and data reduc¬ 
tion, section [3] shows the observational results, which are 
analysed in sectionUin terms of visibilities (section l4.1l) . dif¬ 
ferential phases fsection |4.2II and simple modelling (section 
14.31) . Finally, section [5] includes a discussion, summarizing 
our main conclusions. 


2 OBSERVATIONS AND DATA REDUCTION 

HD 100546 was observed on 2014 April 16 with the 
VLTI/AMBER interferometer. The night was photometric, 
with seeing below 1". Three unit telescopes (UTs) were used, 
covering three different configurations in terms of projected 
baselines (Bs) and position angles (PAs): UT4-UT1 (119- 


123 m, 243-258°), UT3-UT1 (79-83 m, 219-231°), and UT4- 
UT3 (60-61 m, 281-296°). The ranges come from the spread 
in the uv coverage of the data set. The corresponding spa¬ 
tial resolutions, A/2B, are between 2 and 4 mas. The high- 
resolution A-band mode with FINITO fringe tracking was 
used, providing a spectral resolution of ~ 12000 (25 km s“^) 
around the Bry line (AA 2.15 - 2.18 pm). The field of view 
is limited to the Airy disc of each individual aperture, i.e. 
60 mas for the UTs in the K band. Data reduction was 
performed using the AMDLIB-V3.0.8 software provided by 
the Je an-Marie Mariotti Center, following s tandard proce¬ 
dures (ITatulli et al.|[^07al : IChelli et al.ll201ll l. A total of 750 
frames were obtained, from which we took the 20% with the 
best SNR to produce the final averaged observables: fluxes, 
visibilities and differential phases. This provides optimal re¬ 
sults with the smallest noise, but it is noted that the final 
averaged data refers to frames that were taken with a max¬ 
imum time delay of 1 hour, which is reflected by the ranges 
provided above for the B and PA values. The standard star 
HD 101531 (TATll) was observed with the same configura¬ 
tion, before and after the observations of HD 100546, and 
data reduction was carried out in the same way as the science 
target. Instrumental artefacts were corrected for by divid¬ 
ing the observables of HD 100546 by those from HD 101531. 
This procedure also allowed us to remove telluric contami¬ 
nation from the observed fluxes of HD 100546. Two telluric 
(water) lines at both sides of Bry (2.163477 and 2.168687 
pm) were used to perform wavelength calibration. Radial 
velocity correction to local standard of rest was applied. 
In the remainder of this paper all velocities are expressed 
with respect to the sy stemic velocity of 9 km s“^ (LSR; 
iKharchenko et ^1200^ . 


3 OBSERVATIONAL RESULTS 

Figure [T] shows the interferometric observables: Bry line 
fluxes, visibilities and differential phases. The closure phase 
is also measured by AMBER, showing a random noise of ± 
6 ° without any clear trace of signal around Bry. This sug¬ 
gests that the line emitting region is centro-symmetric on 
the sky, within the errors. 

The total Bry flux (top panels of Fig. [TJ shows a single- 
peaked profile with an equivalent width (EW) of -13.6 ± 
0.2 A, where the uncertainties come from the standard de¬ 
viation of the measurements considering the three baselines. 
Because the dusty inner disc extends up to distances < 40 
mas (TATll), the field of view of our AMBER observations 
(60 mas) covers the spatial scale from where all the contin¬ 
uum in K originates. Similarly, the spatial scale from which 
the Bry emission originates is fully included in the field of 
view, as will be argued in this paper. Therefore, the reported 
EW is not affected by possible lacks in the continuum or 
the line fluxes, representing the “real” EW that could be 
measured by any spectrog raph (for a wider discussion, see 
lOudmaiier fc de Wid[2013 1. The line is broadened ~ ± 200 
km s“^. 

Continuum visibilities were fixed to the values provided by 
TATll from K band calibrated AMBER data of HD 100546. 
According to that work (see Fig. 1 in TATll), the visibil¬ 
ity of the continuum in K is close to zero for our three 
baselines due to the (spatially extended) inner dust disc. 
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Figure 1. Br 7 normalized fluxes, calibrated squared visibilities, 
and differential phases of HD 100546 for the baselines indicated. 
Vertical dotted lines show the Br 7 zero-velocity reference point 
for an observer comoving with HD 100546. The profiles have been 
coloured to indicate different spectral channels: 0 ib 20 km s“^ 
(black), 60 zb 40 km (orange), -60 zb 40 km s“^ (light blue), 
130 zb 30 km (red) and -130 zb 30 km (dark blue). Vertical 
error bars are the standard deviations of the adjacent continuum. 
Continuum corrected squared visibilities and phases as derived 
from equations ^ and [2] (see section II are plotted above the ob¬ 
served ones. Their errors were derived from propagation of the 
individual uncertainties. 

The squared visibility increases at wavelengths correspond¬ 
ing to the Br 7 transition (mid-panels of Fig. [T|, indicating 
a more compact emitting region than the continuum. The 
baseline UT3-UT1 (and maybe UT4-UT1, considering the 
errors) shows a double peaked visibility with a central dip 
at 0 km s“^. The UT4-UT3 baseline shows a single-peaked 
squared visibility with a redshifted maximum, although the 
corresponding continunm-corrected visibility (see section a 
is also suggestive of a double peaked profile. 

The phase at Bry wavelengths is different than that of the 
adjacent continuum (bottom-panels of Fig. [T])i indicating a 
spatial displacement of the photocentre of the line emitting 
region with respect to that of the continuum. The differen¬ 
tial phase peaks at different wavelengths depending on the 
baseline considered. 


4 ANALYSIS 

The observed line fluxes, visibilities and differential phases 
(F, V, (f>) include the contribut i on of the continuum (Fc, 
Vc, <(>c). Following I Weigelt et al.l |2oB)i the visibilities and 
phases characterizing the pure line emitting region (Vi,4>i) 
are given by 


and phases are overplotted in Fig.[T] In the following sections 
we will refer to continuum-corrected visibilities and phases. 

4.1 Visibilities: spatial scale of the Bry emitting 
region 

Fig. [2] shows the continuum-corrected squared visibilities 
against the spatial frequencies for the three baselines 
explored. Different velocity bins within the Bry emission 
are represented, using the same colour code as in Fig. [T] 
Our limited uv coverage prevents a detailed modelling of 
the Bry brightness distribution, given that different models 
can reasonably fit the data from only three baselines. How¬ 
ever, we can extract some quantitative information on the 
relative size of the Bry emitting region, compared with that 
obtained from the K band continuum visibilities (TATll, 
plotted with black filled circles in Fig. El). As previously 
mentioned, it is assumed that continuum visibilities are 
mainly representative of the inner dust disc. Therefore, 
the spatial scale of the line emitting region derived below 
is with respect to the inner disc, and not with respect to 
any other source of emission like the central star. This is 
based on two facts: first, the K band flux primarily comes 
from the inner disc, as one can see in the SED shown in 
Fig. 2 of TATll (see also below and sections 14.21 and 14.31) . 
Second, and according to the same work, the visibility of 
the continuum in K is close to 0 for our baselines, again 
indicating that the dominant fraction of the flux comes 
from a resolved source and not from the central star. 
The values in Fig. El can be roughly fitted with simple 
brightness distributions (Gaussian, disc, ring, and elongated 
ring). In all cases the inferred Bry spatial scale is smaller 
than that of the inner dust disc, typically by a factor ~ 
0.8. The corresponding spatial scales associated with the 
red/blueshifted line emitting regions are smaller. It is noted 
that the line visibility level barely decreases from ~ 37 to 
55 Mrad . Interpreting the line visibility as the ratio 
between the correlated, unresolved flux and the observed 
one. Fig. El suggests that at ~ 55 Mrad (spatial scales < 
1.8 mas) there is an unresolved component contributing ~ 
35% to the total line flux. The Bry visibilities of the central 
spectral bins are actually more consistent with a brightness 
distribution that includes a point source contributing by 
that amount, plus a disc contributing ~ 65% to the total 
flux (solid line in Fig. El see sections 14.21 and 14.311 . The disc 
component was assumed to have the same inclination and 
PA as the dust inner disc, and extends up to its inner rim 
(i ~ 33°, PA ~ 140°, ~ 0.25 an; see TATll). The 

best model for the continuum visibilities from TATll is 
also overplotted (dotted line). 


4.2 Differential phases: photocentres of the Bry 
emitting region and the continuum 


F?Vf = -f F'^Vf - 2FV ■ EA4 • cos <(>, 


sin (j)i 


sin 



( 1 ) 

( 2 ) 


with Fi = F - Fc- The values for F;, and Vc are the averaged 
ones at both sides of Bry. Continuum-corrected visibilities 


Differential phases provide information of photocentre dis¬ 
placements along the direction of the baselines, on angular 
scales that can surp ass the nominal resolution of the inter¬ 
ferometer. Following iLe Bouauin et al.l ll2009lf . wavelength- 
dependent photocentre displacements can be derived from 
the continuum-corrected differential phases {(j>i) by 
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Figure 2. Continuum (black solid circles), and continuum- 
corrected Br 7 squared visibilities (open circles; different colours 
representing the Brq spectral bins, as in Fig. (TJ versus the spatial 
frequency. The dotted line represents the fit for the continuum vis¬ 
ibilities. Following TATll, the brightness distribution combines a 
point source contributing 25% to the total flux, and an inclined 
ring contributing ~ 75% {i ~ 33°, PA ~ 140°, r^i^ ~ 0.25 au). 
The solid line represents the fit for the Brq visibilities of the cen¬ 
tral spectral bins. The brightness distribution combines a point 
source contributing 35% to the total flux, and an inclined disc 
(with the same PA and inclination, and limited by the dust inner 
rim) contributing 65%. 



Figure 3. On-sky position (angular distances were transformed 
to au assuming a distance to the system of 100 pc) of the dusty 
inner rim (dashed line, from TATll: PA = 140°, i =33°, Rd = 
0.25 au), and an upper limit for the disc truncation radius limit¬ 
ing a possible MA region (dot-dashed line). The 2D-photocentre 
solutions (from a reference frame centred on the star) of the con¬ 
tinuum are the black symbols inside the disc truncation radius 
(representing 360 velocity bins homogeneously distributed at both 
sides of the Bry line), and those of the Bry emission are the cir¬ 
cles (same colour code as in Fig. [!}. The top left vectors show the 
directions of the baselines. The embedded panel is a zoom-in of 
the innermost regions. 


where pi is the projection on the baseline Bi of the 2D- 
photocentre vector. The 2D-photocentre vectors at different 
wavelengths can therefore be recovered by solving the above 
system of equations. As a hrst order approach (see next sec¬ 
tion), it is assumed that the dust continuum emission is 
homogeneously distributed along the inner disc. Therefore, 
the symmetry of the system forces its photocentre to be 
centred on the (unresolved) star, which is initially assumed 
to be the reference frame from which photocentre displace¬ 
ments are measured. Fig. [3] shows the on-sky positions of 
the derived 2D-photocentre displacements with respect to 
the central star, considering spectral bins that represent the 
continuum and the different components of the Bry line. The 
dusty inner rim is represented by the dashed line (TATll). 
The dot-dashed line is an upper limit for the disc truncation 
radius (Rt ~ 0.01 au), i.e. the maximum stelleocentric dis¬ 
tance from which gas can be accreted magnetospherically. 
This has b een derived from th e maximum magnetic field 
pr ovided inlHub rie et al l (l200fll) . and the expression for Rt 
in iTambovtseva et al.l (12014 ) (see section [5]). 

The 2D photocentre solutions of the continuum bins 
are oriented in the SE-NW direction, which is the same 
as the major axis of the dusty disc reported by TATll 
(PA =140 ± 16 °). A linear fit to our data confirms this 
result, providing PA = 139 ± 6°. In addition, the 2D 


photocentres of the different spectral bins of the line are 
also aligned (PA = 152 ± 8°) with the major axis of the 
system. Given that our uv coverage is not homogeneous, it 
cannot be ruled out that those alignments could result from 
the similar directions of the three baselines covered here. 
Other baselines should be explored to conhrm those results, 
preferably oriented in directions perpendicular to the ones 
analysed. In addition, the 2D photocentre solutions of the 
redshifted parts of the line are located towards the NW 
direction, and the blueshifted parts towards the SE. This 
cannot result from the previously mentioned configuration 
of the baselines, and constitutes an indication of a rotating, 
disc-like structure. Interestingly, the sense of rotation is the 
same as previously reported from other gas tracers in outer 
parts of the disc of HD 1005 46 (i.e. counter-clockwise; see 
lAcke &: van den Anckeil |200 b) . Moreover, there is a rough 
tendency for the high-velocity bins (red and blue circles) to 
be closer to the 2D photocentre of the peak emission (black 
circles), and for the low-velocity bins (orange and light blue 
circles) to be further away. This is suggestive of Keplerian 
rotation, although the error bars are large to provide a 
firm conclusion based on the observations alone (see the 
next section). Finally, the 2D photocentre of the Bry peak 
emission is displaced ~ 0.03 au NE from the continuum 
photocentre, which in this section has been assumed to be 
in the central star. 
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Figure 4. Modelled inner dust rim and velocity surfaces of the 
gaseous disc around HD 100546, as described in the text. North 
is on the top and east is on the left. The NE part is further 
from us and the SW is closer. Symmetric distributions of gas at 
velocity ranges -86 —73, -33 —20, 31 - 44, and 84 - 97 km 
s“^ are represented in the top left, top right, bottom left and 
bottom right panels, respectively. The colour gradient indicates 
the relative brightnesses in each panel, with the darkest colour 
representing the maximum flux. 



-500-400-300-200-100 0 100 200 300 400 500 


v (km s *) 

Figure 5. Observed Bry flux and modulus of the photocentric 
displacements along the three baselines (solid lines with error 
bars) versus model results (dashed lines) as described in the text. 
The modelled flux consists of a Gaussian plus a Keplerian profile 
(dotted and dash-dotted lines in the top panel). 


4.3 Inner dust and gas modelling 


In order to further explore the photocentric shift between 
gas and dust emission we performed simple modelling us¬ 
ing a Keplerian rotating gas plus a dusty inner disc. As 
suggested in section [ATI the Bry flux consists of two com¬ 
ponents, for which we assume a rotating disc and an unre¬ 
solved Gaussian emission region. The rotating component is 
limited by the dust inner disc at 0.25 au (TATll), and an 
assumed i nner hole at Rf. Although a disc wind may play 
a role (e.g. iGarcia Lonez et al.ll2015l ). this is not considered 
in our model. The wind contribution would be resolved out 
and therefore it would not contribute to the visibilities and 
differential phases available so far. A continuum image of the 
in ner dust rim was crea ted using the radiative transfer code 
in I Whitney et all (l2003l ). which has been previously used to 
mo del high angular res olution data of young star discs (see 
e.g. lAkeson et al.ll200a ). The dust structure is a standard, 
flared accretion disc with H{r) = Ho{r/ro)^ {P = 1.25), 
and the inner rim is assumed to be vertical (as opposed to 
curved). The inclination and PA for both the dusty inner 
disc and the Keplerian, Bry emitting region are the ones de¬ 
rived by TATll. 

Fig. H shows the result of this modelling consider¬ 
ing different velocity bins for the gas. From the model 
and the inclination of the source, the inner rim in the 
NE region -further away from us than the SW- dom- 
ina tes the total K band contin uum em ission (see also 
e.g.lAcke fc van den Anckeil ll2006l) . TATll. fAvenhaus et al.l 
ll2014h . Fedele et al.l ( 201^ . which results in a continuum 


photocentric shift of ~ 0.03 au from the central star in that 
direction. The photocentre displacement in each spectral 
channel of the Bry emission was calculated and converted 
to phase shifts by projecting on to the three VLTI base¬ 
lines used. The modelled photocentre displacements and flux 
were convolved with a Gaussian point spread function for a 
spectral resolution of 12000. The results from the model are 
compared to the observations in Fig.[S] The modelled photo¬ 
centre displacements (bottom panels) are in good agreement 
with those observed. The peaks in the photocentre displace¬ 
ments correspond to the peaks of the flux component related 
with the disc’s velocity profile, which also coincide with the 
two peaks shown in the visibilities of Fig. [T] The modelled, 
correlated (unresolved) flux is smaller than the observed one, 
which is consistent with the low observed visibilities. Clearly, 
a significant fraction of the total line flux is resolved out at 
the baseline lengths employed. The modelling also reveals 
that large inner rim heights cause stronger photocentric dis¬ 
placements, whereas other parameters such as the inclina¬ 
tion or the inner radius for the gas have a smaller effect. 
A minimal scaleheight {H/r ~ 0.02) at the inner rim was 
found to best fit the differential phases, consistent with the 
idea that the dusty inner rim is not puffed up (TATll). 

In summary, the continuum photocentre of the star-disc 
system does not coincide with the central star -as it was 
assumed in the previous section- but is displaced 0.03 au 
to the NE. In addition, although the observed Bry profile is 
single peaked, the correlated flux is double peaked and the 
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phases indicate that the resolved 617 emission is in Keple- 
rian rotation and interior to the dusty inner rim. 

The fact that the continuum photocentre is displaced ~ 
0.03 au NE from the star results from the inner disc orienta¬ 
tion, with the furthest part directly illuminated by the star. 
Since the continuum photocentre is the reference frame from 
which the Br 7 one is measured, and given that this is located 
an additional ~ 0.03 au NE from it (previous section), the 
photocentre of the gaseous disc traced by Br 7 is shifted by 
~ 0.06 au from the central star in the same direction. The 
most plausible explanation for this shift is that the gaseous 
disc is flared, implying that the (furthest) NE part of it of¬ 
fers a larger projected surface than the SW. This interpreta¬ 
tion is analogous to the one given bv lAcke fc van den AnckCT 
ll2006ll in order to explain the photocentric shift to the NE 
also observed in the [01] 6300 line. Despite this type of mod¬ 
elling is out of the scope of this work, a simple estimate of 
the flaring angle can be derived from geometrical consider¬ 
ations. A photocentre displacement of ~ 0.06 au to the NE 
requires the projected area of the NE disc to be ~ 2 times 
the one of the SW. For a constant flaring angle a (defined 
as the angle subtended by an inclined disc with respect to 
a flat disc), the ratio of the further-to-nearby projected disc 
surface equals cos {i - a) / cos {i + a). For i ~ 33°, this 
requires cr to be ~ 25°, or H/r — tan(a) ~ 0.50. 

Alternatively, it could be argued that the SW part of 
the gaseous disc is partially obscured by the SW inner rim. 
However, this would need the inner rim to be significantly 
puffed-up, for which there is no evidence. Finally, an exciting 
possibility that could eventually explain the Bry photocen¬ 
tric shift is that the Keplerian gaseous disc is not centred 
on the star but on a different source displaced ~ 0.06 au 
NE from it. This disc could be associated with an accret¬ 
ing, unresolved sub-stellar/planetary companion. The corre¬ 
sponding rotational period around the star would be around 
10 days, for which the probability that we have detected it 
precisely when that is located in the same direction than 
the bright inner rim would only be around 10-15 %. Al¬ 
though this possibility cannot be completely discarded with 
the data in hand, we consider the flared Bry disc centred 
on the star as the most likely scenario. Additional baselines 
and uv-coverage are necessary to carry out a more detailed 
modelling that allows us to confirm/discard this view. 


5 DISCUSSION AND CONCLUSIONS 

We have presented the first high-resolution Br 7 spectro- 
interferometric data of HD 100546. Our results are 
consistent with a Brq emitting region inside the inner rim 
(< 0.25 au) and extending very close to the star. Most 
of the Br 7 emission can be associated with a Keplerian, 
disc-like structure coplanar with the dusty inner disc, and 
rotating counter-clockwise (SE towards us and NW moving 
away from us). Simple modelling is also consistent with this 
view. The Brq photocentric shift of ~ 0.06 au NE from the 
star is most probably due to the flared geometry (~ 25°) 
of the gaseous disc. Additional baselines and uv coverage is 
needed to fully understand the nature of the Bry disc. 


5.1 Accretion in HD 100546 

The Bry line is a common accretion tracer in PMS stars. 
The (non-photospheric) Bry luminosity of HD 100546 can 
be derived su btracting a photospheric EW of ~ 24 A (from 
iKurucj I 1993 I) to the observed EW, and then assuming 
that the disc continuum flux in K represents 75% of the 
total observed flux (TATll). This is in turn given by a 
2MASS typical K magnitude = 5.42. The Brq luminosity 
can then be associated with a mass ac cretion rate from 
the ac cretion-line luminosity relation in [ Mendieuti'a et al.l 
ll201lll (also from iDonehew fc BrittainI 1201 l| j. which was 
derived assuming MA operating in HAeBe stars. The mass 
accretion rate of HD 100546 obtained this way is ~ 3.6 
X 10“^ M© yr“^. The same result within uncertainties 
is obtained from MA mo delling of the UV-Balmer excess 
by iFairlamb et al.l (l2015l l: 2.5 x 10“^ M© yr ^. When 
the UV excess is instead modelled from BL accretion, 
a higher value up to 2 x 10“® M© y r~^ is obtained 
i Bjondel fc Tiin A Diid I 2 OO 6 I I. However, iMuzerolle et al.l 
I 2004h showed that accretion rates larger than 10“® M© 
yr”*^ would make the inner gas optically thick. As discussed 
in that work, that would have a double effect: first, the 
gas would significantly contribute to the total emergent 
flux in the K band; secondly, the stellar radiation would 
not heat the inner dust rim, which would therefore be 
located closer to the star than the dust sublimation radius. 
In contrast, the K band flux in HD 100546 is dominated 
by the inner dust disc, and the location of the inner rim 
is consistent with the dust sublimation radius (TATll). 
Therefore, the mass accretion rate estimated assuming BL 
llBlondel fc Tiin A Diid l2006h is too high to explain the 
observed inner disc properties of HD 100546. In addition, 
the inner rim is not puffed-up according to TATll and 
our modelling. Consequently, the amount of inner gas we 
have detected should be large enough to block the stellar 
radiation and prevent the inner rim to be heated and 
puffed-up, but at the same time small enough to prevent 
the presence of dust at distances closer than the dust 
sublimation radius. 

There are additional lines of evidence suggesting that 
accretion in HD 100546 could be magnetospheri c , and 
not through a BL. Equation 3 in Ijohns-Krull et all (Il999lj 
provides a lower limit to the magnetic field necessary to 
drive MA, in terms of the stellar parameters and accretion 
rate. Assuming a mass accretion rate of ~ 10“^ M© yr“^, 
the minimum stellar mass and maximum stellar radius 
allowe d by the uncertainties provided by IFairlamb et al.l 
ll2015ll . and a minimum rotation period of 0.26 days (from a 
maximum projected rotational velocity vs ini = 65 km s~^ 
and a minimum inclination i = 22°, from iGuimaraes et al.l 
l2006l . and TATll), HD 100546 could require a magnetic 
field of only several tens of Gauss to drive accretion 
magnetospherical l^Q. This is c o nsiste nt with the magnetic 
field measured bv iHubrig et al.l ll2009lj . 8 9 ± 26 G, although 
non detections have also been reported dPonati et al.lll9^ : 

^ Note that large uncertainties are involved in the calculation of 
the magnetic field necessary to drive MA. For the typical error 
bars of the stellar parameters of HD 100546, that ranges between 
a few tens of Gauss and a few kilo-Gauss 
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iHubrig et al.l[201^ . In addition, if accretion in HD 100546 
is actually magnetospheric the Keplerian gaseous disc 
should be truncated at a stelleocentric radius that increases 
with the stellar magnetic field and radius , and decreases 
with the accretion rate an d stellar mass IIElsner fc Lartibl 
Il97lfl . Using equation (6) in iTambovtseva et alj (l2014fl . the 
disc truncation radius of HD 100546 should be < 0.01 au 
(Fig. |3]). This value is consistent with the Keplerian radius 
inferred from the width of the Bry line profile (~ ± 200 
km s“^; section El); once this value is de-projected using 
the inclination in TATll, the corresponding Keplerian 
distance is < 0.02 au. It is noted that, in contrast to 
optically thick lines lik e Ha, Bry is mainly b roadened 
by the Doppler effect dTambovtseva et al.l l2014l l. If the 
Keplerian disc is truncated by the stellar magnetic field the 
gas would then fall ballistically on to the stellar surface. 
This could eventually be traced from the presence of 
redshifted absorptions at free-fall velocities in the profiles 
of several lines. In fact, previous spectroscopic analysis 
involving optic al/UV lines suggest MA/ejection pro cesses 
in HD 100546 llVieira et al.lll9^ : iDeleuil et al.ll2004l ~). with 
redshifted absorptions at velocities comparable to free-fall 
dCuimaraes et aHl2006l ). However, those signatures are not 
observed in the Bry profile of HD 100546. The unresolved 
component contributing ~ 35% to the visibility (section 
14.111 could be related to magnetospheric infall, but the 
small spatial scales involved cannot be probed from our 
observations. 


5.2 Origin of the inner disc gas 

With a mass accretion rate of a few 10“^ M© yr“^, and 
a gaseous inner disc mass of 10~® M© (from the dust 
inner disc mass in TATll and a gas to dust ratio of 100), 
the inner disc would be depleted of gas in le ss than a year. 
Lower values for the gas to dust ratio (e.g. iBrittain et al.l 
I 2 OO 9 I : iBenistv et al.llioiol ') or differences within 1 dex for the 
accretion rate would not significantly modify that result, 
which indicates that the survival time of the inner disc gas 
can only be a few months/years. This would require that the 
observed gas in the inner disc is being replenished, perhaps 
by flows from the outer to the inner disc. This type of 
planet-boosted fl ows has been observe d in similar HAeBes 
like HD 142527 (ICasassus et al.ll2013ll. which is accreting 
at a similar rate than HD 100546 ( Mendieutfa et al. I [ 2 ^. 
If planets are indeed being formed in this type of systems, 
they do not seem to have a significant effect on the stellar 
accretion rates, gi ven that 10~^ M© yr~^ is a typical value 
for HAeBe stars (iMendiguti'a et aDl201lf) . In other words, 
despite planets in formation ca n induce gas transfers be¬ 
tween different parts of the disc llDodson-Robinson fc Saiv^ 
[2 oT 3), the amount of gas trapped by them would be 
low enough to keep the stellar accretion rate practi¬ 
cally unaltered, as expecte d from the compa r atively low 
plane tary accretion rates llClose et al.l I20l4 IZhou et al] 
I 2 OI 4 IL It is also noted that the fact that the accretion 
rate associated to the Bry emission is several orders of 
magnitude larger than for sub-stellar/planetary objects 
argues against the unlikely explanation that the ~ 0.06 au 
Bry photocentric displacement could be associated to an ac¬ 


creting disc centred on such a close companion ('section l4.3l) . 


5.3 Final remarks 

As argued before, and regardless from a possible MA con¬ 
tribution to the Bry line, most of its em ission comes from a 
more extended, Keplerian rotating disc. iKraus et al.l ll2008l j 
found a correlation between the Ha profile and the extent of 
the Bry emitting region with respect to that of the contin¬ 
uum characterizing the inner dust rim, which in turn would 
be related with the physical origin of the line. Stars with a 
P-Cygni Ha profile show particularly compact Bry-emitting 
regions {Kbt-i/^ cont < 0.2) more consistent with accretion, 
while stars with a double-peaked or single-peaked Ha pro¬ 
file show a significantly more extended Bry-emitting region 
(0.6 < Rflr 7 /Rcont < 1.4). The Ha profile of HD 100 546 
is variable, but always double-peake d (I Vieira et al.ll 1999 1 or 
single-peaked (iFairlamb et al.n2015l ). This, along with the 
relative size of the Bry emitting region /^cont ~ 0.8; 

section [441 . again suggests that the major contribution to 
the Bry line in HD 100546 is not directly tracing accre¬ 
tion. Moreover, MA modelling of the Bry line considering 
stellar parameters similar to those of HD 100546 produces 
broader, double-p eaked profiles that are no t consistent with 
our observations (iTambovtseva et al.ll201'i l. We emphasize 
that this type of modelling assumes a disc truncation ra¬ 
dius much larger than the one derived for HD 100546. This 
can strongly affect the broadening of the modelled line, as 
well as its double/single-peaked nature, depending on the 
spectral resolution assumed. Smaller magnetospheres could 
indee d be typical for HAeBe stars (ICaulev fc Johns-Krulll 
l2014h . Generally speaking, previous models that reproduce 
the Bry emission exc l usive l y from magnetically channelled 
gas (iMuzerolle et al.l 1200 ll : [Kurosawa &: Romanova! l2013f ) 
are not consistent with our observations. 

In summary, whereas accretion in HD 100546 could be mag¬ 
netospheric, the major contribution to the Bry emission 
comes from a Keplerian, probably flared, gaseous disc, and 
not from magnetically channelled gas. The accretion rate 
and the inner disc mass indicate that the observed gas should 
be replenished, perhaps from planet-induced flows from the 
outer to the inner disc as have been reported in similar ob¬ 
jects. The peculiar characteristics of HD 100546 make it the 
perfect laboratory both to test planet formation and to un¬ 
derstand how accretion proceeds in HAeBe stars. 
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